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Abstract 

The matrix 2 x 2 spectral differential equation of the second order is considered on x in (—00, +00). We 
establish elementary Darboux transformations covariance of the problem and analyze its combinations. 
We select a second covariant equation to form Lax pair of a coupled KdV-MKdV system. The 
sequence of the elementary Darboux transformations of the zero-potential seed produce two-parameter 
solution for the coupled KdV-MKdV system with reductions. We show effects of parameters on the 
resulting solutions (reality, singularity). A numerical method for general coupled KdV-MKdV system 
is introduced. The method is based on a difference scheme for Cauchy problems for arbitrary number 
of equations with constants coefficients. We analyze stability and prove the convergence of the scheme 
which is also tested by numerical simulation of the explicit solutions. 


Introduction 


lere are two complementary approaches to integrable systems: analytical and numerical ones to be bevel¬ 
ed. Even most profound analytical 1ST method cannot give explicit solution of general Cauchy problem 
while numeric can, but is rather compulsory in use: calculations could need powerful computers. May be 
most transparent of analytical methods are based on algebraic structures associated with a problem. To 
such structures belongs Darboux transformations covariance of Lax representation of nonlinear equations 
that yields a powerful tool for explicit solutions production. We investigate applications of special kind of 
such discrete symmetry - to be called elementary ones m- Its elementarity simply means that a product 
of such transformations generate the standard one 0,|T1. Here we study combinations of such transforms 


that do not coincide with binary ones |]^ and hence are not so known. 

The main ideas of numerical integration of such integrable systems go up to the famous properties of 
the equations as the Lax pair and inhnite series of conservation laws existence [^. From a point of view 
of general theory of such systems some hopes are concerned with a development of the hnite-difference or 
other approximations of the systems. Namely if one could prove a convergence and stability theorems for 
such difference systems (existence of solutions is implied), a way to existence and uniqueness of solutions 
is opened |P. 

The coupled KdV-MKdV system arises in many problems of mathematical physics. Some integrable 
systems are associated with a polynomial spectral problem and have Virasoro symmetry algebras are 
considered [^. A dispersive system describing a vector multiplet interacting with the KdV held is a 
member of a bi-Hamiltonian integrable hierarchy [§]. Recently a multisymplectic numerical twelve points 
scheme was produced. This scheme is equivalent to the multisymplectic Preissmann scheme and is applied 







to solitary waves over long time interval [1^]. The conpled KdV-MKdV system is also connected to other 
physical applications O . 


The general system we consider in this work have the following form 

+ E = o, 

m,k 

n,l 


( 1 ) 


where n, m, k =1,2,...N are the dependent variables nnmbers. Nonlinear coefficients are I =1,2,..5 
and dn are dispersion coefficients. 

In particular, for the system under consideration (iV = 3) the variables 0^,9"^, 9^ are denoted by f, u, v 
to have the integrable system 


10 


ft + 1/... +1(«/). - luf = 0 

_3 3 

^^xxx 2 ^^^ -|- ^^xj ^yjx'^jx 0 

1 3 3 3 3 

T f^'^xxx T 2^9;^ ~^^xxf T ~Uxfx 0 


( 2 ) 


The Lax pair is given in [|T^]. The system exhibits two integrable reductions having explicit solutions. 


Hirota-Satsuma ||, and a two components KdV-MKdV system |^, |^. Krishnan ^ showed that a 


generalized KdV-MKdV system have solitary wave solutions and investigate the effects of increasing the 
nonlinearity of one variable on the existence of solitary waves. Some form of KdV-MKdV system have 
explicit solutions in terms of Jacobi elliptic functions 0]. 

In this work we present explicit solutions for a system of three equations (^) that have not been specihed 
TU|, O. We study this two-parameter explicit solutions and show effects of choosing these parameters 


m 


on the solutions. We demonstrate the use of two arbitrary elementary DTs [^| and its special choice that 
holds a hereditary of the reduction to built explicit solutions to the KdV-MKdV system (^. 

Also we modify a numerical method ||^, for solution of system (|l|). It is a difference scheme for 
Cauchy problems for arbitrary number of equations with constants coefficients. The scheme preserves two 
conservation laws for the KdV type equations and the order of error of the difference formulas is improved 
0, H- The convergence is proved and stability is analyzed giving the conditions taken in account in 


choosing time and space step sizes 

The present work is organized as follows. Section 2 introduces the matrix spectral equation of the 
second order with 2x2 matrix coefficients and two elementary DTs. We select the second equation of 
the Lax pair and derive the compatibility conditions. The product of these two transformations yields 
the standard DT |^, |T^. Section 3 illustrates how the, hrst, elementary DT is used to produce solution 


to the KdV equation as well as the general evolution equation generated by the compatibility conditions 
of the Lax pair. Explicit solutions are introduced for the case of zero initial potentials of the matrix 
problem. In section 4 we consider a reduction constraints on the potential of the matrix spectral equation. 
This reduction gives an automorphism that relates two pairs of solution of the spectral equation for two 
spectral parameters. We use this results in the compound elementary DTs to produce an explicit solution 
to a coupled KdV-MKdV system that results from the compatibility conditions of Lax pair under this 
reduction. The effects of these parameters on the solution (reality, singularity) is analyzed. Section 5 
introduces a numerical method for solving coupled KdV-MKdV system (|). We produce a difference 
scheme for a Cauchy problem with initial condition rapidly decreasing at both inhnities. The main steps 
of the scheme convergence and stability analysis is shown while the details are explained in appendix A 
and B. The scheme is tested by applying it to integrable coupled KdV-MKdV system and the numerical 
results are compared with explicit formulas obtained in section 4. 




















































2 Lax pair spectral equations and the elementary DTs 

Consider a matrix spectral equation of the second order with spectral parameter A and 2 x 2 matrix 
coefficients. 

+ + (3) 

where the vector d' = (d'l, \l' 2 )^and the matrix potentials are U = {wp}, F = {fij, fa = 0}, i = 1, 2 while 
(73 = diag{l, —1) is the Pauli matrix. 


For equation @ we perform two elementary Darboux transforms 


d'l — + eiid'i + ei2d'2 , eii — — ( ^pix + -fi2'^2 



The hrst one is 
1, ei2 = /12/2, 


'I /2 — T 2 + 621 d'l , ^21 — — </32/</9i, 


+ eil) + ei2 </?4 , ^2 = + €21 </ 33 - ( 4 ) 

where {ipi, (^ 2 )^ and ((^ 3 , (^ 4 )^ are two solutions of (H) corresponding to different spectral parameters. 

Substituting the above expressions for Ti, ^2 info (|^) and collecting the coefficients of Ti, \['2 and their 
derivatives we obtain the expressions for new potentials as 


f 12 — ni 2 + fi 2 € 11 , 


/21 — “2 621 , 


Wii — Uii — 2eiix — f 12 ^21 — /21 € 12 , 

ni2 = Ui2x — ^12xx + €iiMi 2 — €l2 + U 22 ) , 

■^21 = / 2 I — 2 e 21 a: — / 2 l€ll, 


U22 — U22 — €21^12 — ^21612 — /2i€12x- ( 5 ) 

The second elementary DT is performed after the hrst one and can be obtained by reversing the indices 
1 —>■ 2 and 2—1 to get, for example, the following potentials 


/2I — “21 + f 21^22, €22 — — 1 ^ 9 ^ 23 ; + 2'^2iV^iJ /V ^25 

rvj /*sj ^ rsj r\_/ . 

U22 = U22 — 2e22x — 121^12 — f 12^21, €12 = — €21 = /2l/2. 

The spectral equation (||) is considered as the hrst equation of the Lax pair, take the second as 

= ^xxx + + CT 

where B = ^diagU + |Fj. + 

and C = |t 4 — ^diaglJx — |(/i2“2i + /2i“i2)-f + |(/i2,x/2i ~ fi 2 f 2 i,x )(^3 + |(“ii ~ U22)cy^F. 


( 6 ) 

(7) 


Equation (|^ is also covariant under transformations (^, (^. The compatibility conditions have the 
following form 

Ft — F^x + B 2 X — 3 U 2 X + 2C*a; + FBx — (j^Ba^Fx -\-UB 
—azBa^U + FC — a^Ca^F = 0, 

Ut — Usx + C 2 X + UC — a^Ca^U + FCx — (r^Ba^Ux = 0. 

and the transformations (^, determine a discrete symmetry of (|^) 


( 8 ) 




3 Solution of two coupled KdV-MKdV equations and KdV 
equation via the first elementary DT 

For a spectral parameter A and a seed potential F, U we obtain the solntions to the pair i). ©■ 

Then performing the first elementary DT to obtain the new potentials F, U which are solntions to the 
system (H). For the case of zero seed potential the solntions, tpi and (p 2 of the system 0.0 have the 
form 


(^2 = (9) 


where Ci, C 2 , di, d 2 are arbitrary constants, a = 
the only nonzero elements ) to the following 

f21t + —f21xxx 


^/X and i is the imaginary nnit. System (j^) rednced (for 


+ -^f2lUUx 


3 

--U 11 U 21 , 


_1 _n 

1 3 3 3 ^32^ 

'^21t + -^U21xxx + ^'^21^113; + -U2lxUll — ^ Ullj21xx + ^'^11/21- 


( 10 ) 


where /12 = 0, M 12 = 0, U 22 = 0 and tildes are omitted for simplicity. This system with explicit solntion 
obtained from (^), (j|) as 

/21 = (^d2e2“'‘ + / (c2 + , 

nil = (8a2ciC2e2“(“'*+")) / (c2 + , 

U2i = (^d2e'“'' - die'^")) / (^C2 + Cie2“(“'*+")) . (11) 

where ci, C 2 , di, d 2 are arbitrary constants and a = a/A. 


The second eqnation in (^) is the KdV eqnation while the remaining are a two components conpled 
KdV-MKdV system that was solved by elementary DT. 


4 Solution of three coupled KdV-MKdV equations via the com¬ 
pound elementary DTs 


Existence of different kinds of antomorphism causes special constraints 
to have 

+ O-lFTa, + CJit/T = ActictsT (12) 

but (JiCTa = —a^ai and consider the conditions aiF = Fai and aiU = Uai that means 



10]. Multiplying (g) by cxi = 


fl2 — /2I — /, Uii — U22 — U, U12 — U21 — V. 


( 13 ) 


So (|T^ becomes 


^ ^ (^1^) = “-^^3 (c^l^) 






The above automorphism T (A) <— criT (—A) relates two pair of solutions ((^ 1 , 992 ) and ( 993 , (^ 4 ) of (j^) 
corresponding to different values of spectral parameter A, —A as 


V93(-A) 

(P4(-A) 


= CTl 


^ M(A)\ 
\(P2(A)y V¥^i(A)y 


Using this result in the elementary DTs (^, (|^) and (|^) to obtain the expressions for the new potentials 
f, u, V. In the case of zero initial potentials these new potentials have the following forms 


u = 


f = 2 


9 ^i( 9 ^ 2 )x - 9 ^ 2 ( 9 ^i)x 


- (^ 2)2 


((Pl)2 - ((P 2 ', 

"" V (v^i)^ - (V^2)2 


+ 2 


/ (flMx - ^2{^l)a 




{^l{^2)x - ^2{^l)x) {{^l)x - {Yl)x) 


+ 


(((pi)2 - 


(14) 


where 994 , 992 are as in (H) with Ci, C 2 , di, ^2 are arbitrary constants and a = y/X. 

The above expressions are solutions of system (|) that reduced under the reduction conditions (|T^ to 
system (|^). 

The choice of the arbitrary constants (ci, C 2 , di, d 2 ) affects on the behavior of the solution in formula 
(|). For example choosing equal constants Ci = C 2 = di = d 2 = 0.5 (we choose the value to be 0.5 to 
simplify the resulting formula but the idea valid for any value) the solutions have the form 

/ = 2 a{sinrjicoshrj 2 — cosriisinhri 2 ) / {cosh‘s ri 2 — cos'^rji), 
u = 2a^{sinriiCOshri2 + cosriisinhri2Y / {cosh‘^ri2 — cos^rjiY, 

V = 20^ {cos3T]icoshr]2 — 2sinriisinhr]2{cos2rii + cosh2Tj2 + 2) — cosr]icosh3r]2) 

/ {cosh‘^r]2 — cos'^riiY. 

where rji = a^t — ax, 772 = o,^t + ax, a = y/X is real. 


(15) 


We see that the above expression (|^ is singular at 772 = 0,771 = mi,n = 0 , 1 , 2 ,.... Hence we have 
singularity at {x = ^,t = ^). 

To obtain continuous solutions we can choose ci = C2,di = d2 = r.ci, r is real constant. We again 
choose Cl = 0.5 following the previous concept. So (^ have the form 

/ = 2 ar (cosh 772 sin 771 — cos 771 sinh 772) / (cosh^ 772 — cos^ 771) , 

77 = (1 — cos 2771 + cosh 2772 + sin 2 r]i sinh 2772) / (cosh^ 772 — cos^ 771) ^ , 

V = 2 a^r (((—7 + 6r^ + 2 r^ cos 2771) cos771 cosh 772 — cos771 cosh 8772) — 2(1 + Y 
+r^ cos 2771 + cosh 2772) sin771 sinh 772))/(—1 + cos 2771 — cosh 2772)^. ( 16 ) 


where 0 , 771,772 as in (|T^ 

Choosing this parameter (r) to be r < 1 gives real nonsingular solutions. The above formula (^ is 
built from elliptic and periodic functions so it does not preserve its symmetry but its localized as shown 
in hgures (l.a,b) below. 











Fig.(l.a) Non-singular solutions, f, u and v (r= 0 . 5 ) , a= 2 , t= 0 . 





Fig.(l.b) Propagation of solutions, f, u and v (r= 0 . 5 ) , a= 2 , t=l. 

Fig.(l) The solutions in does not preserve its symmetry but its localized. 

Choosing the parameter r to be r > 1 in formula gives singular solutions as shown figure ( 2 ) below. 



Fig. ( 2 ) Singular solutions, f, u and v (r= 2 ), a= 2 , t= 0 . 

Moreover the choice of these arbitrary constants (ci, C2, di, ^2 ) as well as the spectral parameter A affects 
the reality of the resulting solution. As example for A = m is real and choosing ci = C2 = di = 

^2 = 0 . 5 , we get real solution 

/ = m{cos2(isinh(2 — sinh(2 — sin(^iCosh2C,2 + sinC,i) / {Q .2bcosh2C,2 — 0 . 25 ) (1 — cos2(i) 

where Ci = 2mx + Am^t, (2 = 2mx — Am^t, while choosing ci = C2 = 1 , di = ^2 = 2 give the following 
complex solution 

f = m* {—S{—5sinh(2Cos2(i + 5sinh(2 + 5sin(icosh2(2 — 5 sin(i) — Si{6sinh(2 
+Ssin2(iCosh(2 + 3cos(isinh2(2 + 6 sin(i))/{ 17 cosh 2(2 + 10 -|- 36 cosCiCosh .^2 

—8cos2(icosh2(2 + 17 cos 2 ^i) 


5 The numerical method 


5.1 The difference scheme 

For the coupled KdV-MKdV system (|^) we introduce a numerical (hnite difference) method of solution 
0 , 1 ^ . This scheme is valid for arbitrary number of equations with constants coefficients and of the form 


Qn,j+l _ 071,j 


nk,j _ nk,j nk,j _ nk,j 

^ - + 9 r,l,ki ^7 n ^ - 


r . 2h ‘ ' 2h 

m,k 

nm,j _ nm,j nk,j _ nk,j nk,j _ rynk,j . nk,j 

-a—+ - Yh - 


0k,j _ 0k,j 

I ^n,5 am,j ak,j"i+l '^7-] 

r,L 

ere i and j are the discrete 
dial step. 


'777,k-7 -7 2h ^ -2h^- = ° 

space and time respectively. The time step is denoted by r while 


(17) 





























































5.2 Stability analysis of the scheme 


We prove stability with respect to small perturbations of initial conditions | 16 , It is the boundness of 


the discrete solution with respect to small perturbation of the initial data. We give here the main steps 
while the details are presented in Appendix A. We can write 


dQn,j+l ^ dQn,j 


rpn,j+l rjinj — l 


n <u 


n,o 

r 


where is perturbations of the discrete solution, d 6 ^'° small perturbation of the initial data and 

is a differentiable operator. Stability required the boundedness of 11 i.e ||T^|| is bounded. 

We found that 




a(r, h) = 2 max 

l,n,m,k 


n^l 

9rrL,k 


I I I 
max I \t 7 ^ r \ 

i,m,k 


^ x.i 


+ r[ max 

l,n^m^k 


n^l 

9 m,k 


max ( 1 6 ^ ■ 

i,m,k 


x.i 


^ x,i 


1 

— max 

h l,n,m,k 


n,l 

9m,k 


{ l/)^J 
max I \ 9 - 

i,m,k 




+ —max|d„|]^ 

n 


( 18 ) 


The scheme is stable if a(r, h) < const. We have here a conditional stability. That is we require that 
r —>■ 0 more faster than h —>■ 0 . Namely we need 


T < (constant) . 


5.3 Convergence proof for the scheme 

We prove that the solution of (|^ converges to the solution of (P if the exact solution is continuously 
differentiable one 0 i- We introduce here the main points for the scheme convergence and give the 
details in Appendix B. 

9 l is the difference solution of (|T 7 |) , u^ is the exact solution. Hence the error v^ is given by vj = 6 l—u\ . 
Introducing L2 norm dehned by ||W || = ( XX ^ 

\ i n 

The scheme converges when the norm of that error ||W || ^ 0 as (r, h —> 0 ) 

We found that ||W^^|| < P{M) ( 9 (r + h^), where P{M) is a polynomial in the bounded constant 
^ = 'T' ^e'rXi cind a as in (|T^). Hence the convergence proved. 


5.4 Numerical calculations and test 

The coupled KdV-MKdV system (|^) is solved numerically using scheme (pT|) with initial condition from 
( 1 ) at f = 0 and the results are compared with the explicit formulas (P. The percentage errors are shown 
in the following plots. 
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Fig .3 percentage errors of the numerical solutions relative to the explicit solutions. 















































The results of the test conhrms the validity of the numerical scheme we propose. It also illustrates the 
errors of evaluation that could be estimated by the resulting inequalities of the scheme convergence proof. 


6 Conclusion 


Darboux transformations covariance of Lax representation of nonlinear equations is a powerful tool for 
explicit solutions production. Here we investigate applications of special kind of such discrete symmetry - 
to be called elementary ones. We use these elementary DT to produce explicit solutions for coupled KdV- 
MKdV system. The iteration of DT can be formulated in form of determinants representations |^, . 

A numerical method for general coupled KdV-MKdV system is introduced. It is a difference scheme for 
Cauchy problems for arbitrary number of equations with constants coefficients. We analyze stability and 
prove the convergence of the scheme. The scheme keeps two conservation laws chosen in analogy with KdV 
type equations. Analyzing stability and proving the convergence beside comparing the numerical results 
with explicit formulas allow us to use the numerical scheme to systems with arbitrary coefficients that is 
presumably non-integrable. Obviously the coupled KdV systems are successfully treated by our scheme |^. 

Acknowledgment We thank S. B. Kshevetskii for useful discussions about numerical scheme for 
the problem under consideration. 


References 

[1] Dodd R and Fordy A, On the integrability of a system of coupled KdV equations, Phys. Lett. A89 (1982), 
168-170. 

[2] Guha-Roy C, On explicit solutions of a coupled KdV-MKdV equation. International Journal of Modern 
physics B3(6) (1989), 871-875. 

[3] Halim A A, Kshevetskii S P and Leble S B, On numerical integration of coupled Korteweg-de Vries System, 
Applied Mathematics Letters (to be published). 

[4] Hirota R and Satsuma J, Soliton solution of coupled KdV equations, Phys. Lett. A85 (1981), 407-409. 

[5] Krishnan EV, Remarkes on a system of coupled nonlinear wave equations. Journal of Mathematical physics 
31(5) (1990), 1155-1156. 

[6] Krzhivitski A and Ladyzhenskaya O A, A grid method for the Navier-Stokes equations, Dokl. Akad. Nauk 
SSSR 167 309-311 (Russian); translated as Soviet Physics Dokl. 11 (1966), 212-213. 

[7] Kshevetskii S P, Analytical and numerical investigation of nonlinear internal gravity waves. Nonlinear Pro¬ 
cesses in Geophysics 8 (2001), 37-53. 

[8] Kupershmidt B A, A coupled Korteweg-de Vries equation with dispersion, J. Phys. A: Math. Gen. 18 (1985), 
1571-1573. 

[9] Lankaster P, Theory of matrices. Academic Press, New York, 1969. 

[10] Leble S B and Ustinov N V, Darboux transforms, deep reduction and solitons, J. Phys. A: Math. Gen. (1993), 
5007-5016. 

[11] Leble S B, Nonlinear waves in waveguides. Springer-Verlag Berlin, Germany (1991). 

[12] Leble S B and Ustinov N V, Korteweg-de Vries-modihed Korteweg-de Vries systems and Darboux transforms 
in U-l and 2-|-l dimensions, J. Math. Phys. 34(4) (1992), 1421-1428. 




[13] Matveev V B and Salle M A, Darboux transforms and solitons, Springer-Verlag, Berlin (1991). 


[14] Ping Fu Zhao and Meng Zhao Qin, Multisymplectic geometry and Multisymplectic Preissmann scheme for 
the KdV equation, J. Phys. A: Math. Gen., 33(18) (2000), 3613-3626. 

[15] Shaohong Zhu, A difference scheme for the coupled KdV equation. Communication in Nonlinear Science &: 
Numerical Simulation, 4(1) (1999), 60-63. 

[16] Tannehill J C, Anderson D A and Pletcher R H, Computational fluid mechanics and heat transfer, Taylor &: 
Francis, Washington, (1997). 

[17] Wen-Xiu Ma and Zi-Xiang Zhou, Coupled integrable systems associated with a polynomial spectral problem 
and their Virasoro symmetry algebras, Progress of Theoretical Physics 96(2) (1996), 449-457. 


7 Appendix A 

Stability analysis of the scheme 


We prove stability with respect to small perturbations (because we consider nonlinear equations) of initial 
conditions. Strictly speaking it is the boundness of the discrete solution in terms of small perturbation of 
the initial data. Consider the differential 

190 ^. 1 }, = { 9^4 e-.i 

and dehne the norm \\d 6 ^\\ = ( h 


We can write = Tl 


r n 

doY = Tiv de^-^ = n [t^Y dO'^ 

where is perturbations of the discrete solution, d 6 '^'° small perturbation of the initial data. 

Stability required the boundedness of fl {T^r) i-® 11 ^^II bounded. We calculate T from (^) as follow 


po 


- 5 .- 1 ,.) + Si, (9hi - 


m,k 


2 h L 


n,2 

9m,k 

'Yh 


(e.”4 (5.+.,. - Si-i,) + pa - ea 


11,3 

ym,k 


2 h 


(Cl - C’O (Cl,. - ^*-1,.) + (Cl,. - ^.-1,.) (Cl - Cl 
(^i+l,. ~ + <^i,. (C^l ~ -t- C^l 


n ,5 

9 m,k 

'Yh L' 


n ,4 

9 m,k 

^Yh L 

oT’Di-’ (5.+1.. - 5.-1..)+ 6T-’ ( sa - sa) Si, +9a (9a - 9 a) 5 .,.) 

rd. 


2 h^ 


25j_|_i -|- 2Si_i i~ 2hj_2^j'] 


Rewriting (| 1 ^) in terms of identity (E), symmetric (S) and anti-symmetric (A) matrices 


II 5 '-^'+^II < r max 

l,n,m,k 

||A^+i < Lmax 

l,n,m,k 


n,l 

9 m,k 


n,l 

9m,k 


max ( \ 9 ^f 

i,m,k \ ’ 

max ( 1 9 Y 

i,m,k 


9k,j 

x,i 

9k,j 


-a 4 max |d„|, 


e-’ , 1 -e^. 


where j '^4 = I, 2 , ..A^, Z = l, 2 ,... 5 . 


(19) 


























||TJ+i ||2 = ||(T^'+ 1 )* TJ+ 1 || = i|(E - + S^+^) {E + + ^^'+ 1)11 

<1 + 2 ii^^'+iii + di/ii+iii + ii^j'+ilD^ 

< ^a{T,h)T ^ 


a{T, h) = 2 max 

l,n,m,k 


n^l 

9 m,k 


-— max 

h l,n,m,k 


max ( \ 6 '^f 

i,m,k 
n,l 

9 m,k 




+ t[ max 

l,n,m,k 


n,l 

9m,k 


I I n'^ii 

max 


i,m,k 


x,i 




{ I J I 
max ( \ 9 ^ I 

i,m,k 


9 


k,j 


+ —max \dn\y 

n 


We have here a conditional stability. That is we require that r 
need 


0 more faster than h — 0 . Namely we 


r < (constant) . h® 


8 Appendix B 

The scheme convergence 

We prove the convergence by proving that the norm of the error (between the difference solution and the 
exact solution) vanishes as the mesh is refined. Let 6 l the difference solution of ( 0 ), u^ the exact solution. 
The error vl is given by vl = Of — uj. 

The scheme converges when the norm ||W || —> 0 as (r, h — 0 ) where the norm is defined as ||W|| = 

(ee 

\ t n J 

substitute in P) by 9 l = vl + ul keeping in mind that for 9 l equation (|T 7 |) is O (r + h^) and using 
the operator T defined by 






k,j k,j k,j k,j k,j k,j 

m,j^i+l '^i—1 _|_ 1+—1 


m,k 


2h 


2h 


I ( ( T 

+ 5 'm,fc(W 


2 h 






m,j k,j k,j m,j m,j k,j k,j 


+(n 


mj\2 I n,3 ^'^i+l ‘^i+l '^i-1 . ‘^i+1 ^i+1 ^i-1 


2h 


- \ I / 

■) + 9m,k\ 


2h 


+ 


2h 


2h 2h 

k,j o k,j I k,j k,j o k,j , k,j 

/ m,j ‘^i+l ~ + +_i mj '^i+1 ~ + '^i-1 \ I n n,5 / m,j k,j 

'^9m,k\^i 2}i Oh ' ^ ^9m,k\^i + 


2h 


o I o n,j n,j 

^i+i - ^i-i ^ ^ ^ Ji+2 - 2Vi+i + 2+Ti - Vil2 i 


2h 


2h 


2 h 3 




ir 


So we obtain 


V, 


n,j+l _ 


/ j ir ' ' J'I 


m,k,i 


(20) 


k.j k,i „ k.j k,i „ m.i mA k.j k.j 

Wiicie k i — 2^ym,k^i 2h ^ ym,k\^i ) 2h ^ 9m,k 2h 2h 
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Using Schwartz inequality so p 0 | ) becomes 

||yi+i||<||T^+i|| nil'll+^||/i|| 

< ||T^+i|| IIT^'II ||U^“^|| +r(T^+^ 11 /^'-^ 


< e 


+ r (i 

< e“^^' ||U°|| +M 


,aT(j-l) 


n,l 

9 m,k 


+ e' 


aT{j-2) II 


+ ... + 


||U^+i|r + M 0 (r + h 2 ),M = rf^' 


Using ||U°|| = 0 , the above inequality has the solution 
||y/+i|| < p Q p ^ P(^M) is a polynomial in M. 

Since M is bounded then ||U-^’''^|| —0 as r, h —> 0 and the convergence proved. 













